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I.  INTRODUCTION 


A.  MID-ATLANTIC  BIGHT  OBJECTIVES 

Better  knowledge  of  operating  environments  may  allow  any  Navy  to  improve  the 
use  of  its  USW  systems  in  different  tactical  scenarios.  Even  though  undersea  sound  and  its 
practical  uses  have  been  continuously  explored  since  World  War  II,  testing  of  other 
techniques  are  necessary  in  order  to  obtain  important  sea  environment  information  which 
may  help  in  the  development  of  accurate  Naval  USW  systems.  To  face  this  challenge,  the 
Mid-Atlantic  Bight  (MAB)  experiment  was  conducted.  The  MAB  experiment  was  an 
intensive  field  study  of  the  shelfbreak  south  of  Long  Island  during  the  strongly  stratified 
summer  season  whose  basic  objectives  were: 

1 .  To  provide  detailed  time  series  of  the  oceanographic  conditions  along  typical 
major  acoustic  paths  by  using  moored  acoustical  arrays. 

2.  To  provide  a  diversity  of  frequencies  and  source  locations  from  broadband  SUS 
sources. 

3.  To  determine  the  temperature  structure  of  the  ocean  environment  by  using 
acoustic  tomography. 

This  study  may  provide  the  Navy  an  alternative  and  simple  means  of  determining 
characteristics  in  unknown  or  poorly  known  operating  environments.  (Pickart  et  al.,  1996) 

B.  STUDY  OBJECTIVES 

The  objectives  for  reaching  the  fundamental  goal  of  this  study  are  as  follows: 

1.  To  determine  the  bottom  sound  speed  by  measuring  the  coherence  of  two 
vertically  displaced  sensors  in  the  shallow  water  channel  and  comparing  this 
function  with  one  obtained  from  a  simple  model  equation  over  a  range  of 
frequencies. 

2.  To  determine  the  bottom  density  by  comparing  mode  travel  times  of  computed 
spectrograms  for  a  phone  from  SUS  explosive  data  and  those  obtained  by 
employing  the  concept  of  an  effective  depth  for  an  ideal  waveguide. 

1 


3.    To  determine  the  bottom  attenuation  by  finding  mode  attenuation  as  a  function 
of  range  extracted  from  spectrograms  and  relating  it  to  loss  per  bottom  bounce. 


II.  EXPERIMENTAL/ENVIRONMENTAL  DESCRIPTION 


The  summer  Mid-Atlantic  Bight  experiment  consisted  of  a  cruise  from  1 9  July  to  9 
August  1996  aboard  the  R/V  Endeavor  from  the  University  of  Rhode  Island.  The  cruise 
included  three  major  components:  acoustic  tomography,  SEASOAR  operations,  and 
hydrography/tracers.  An  AXBT  flight  was  also  conducted  during  this  time  period,  as  well 
as  a  SUS  drop  fly-over.  The  acoustic  portion  of  the  experiment  consisted  of  three 
components: 

1 .  The  moored  acoustic  transmission  array,  or  tomography  array 

2.  The  moored  physical  oceanography  array,  which  provided  detailed  time  series 
of  the  oceanographic  conditions  along  our  major  acoustic  path 

3.  The  broadband,  explosive  (MK61  SUS  shot)  source  component. 

Multiple  lines  of  SUS  drops  were  deployed  in  the  area  by  P-3  aircraft  and  monitored  by  two 
vertical  line  arrays  (VLAs)  moored  on  the  shelf  in  ~85m  water.  (Pickart  et  al.,  1996) 

To  achieve  the  objectives  of  this  thesis,  eleven  broadband  SUS  shot  transmissions 
recorded  at  the  NW  VLA  were  selected.  All  of  these  were  dropped  along  a  single  track 
from  a  point  at  latitude  40°16.2'N,  longitude  71°10.2'W  to  a  point  at  latitude  40°16.8'N, 
longitude  70°57.0'W.  The  sampling  rate  was  3906.25  Hz,  the  duration  of  one  record  was 
8.3886  sec  and  the  number  of  samples/record  was  32768.  The  arrangement  and  the  location 
of  the  acoustical  elements  used  for  this  study  are  shown  in  Figure  1 .  This  region  of  the 
experimental  site  has  a  nominal  water  depth  of  ~  90  m  with  little  variation  (<10m)  from  the 
VLA  to  the  positions  of  the  selected  SUS  drops.  The  receiving  VLA  had  16  hydrophones 
of  which  only  the  upper  8  successfully  recorded  the  SUS  transmissions  with  an  upper  depth 
of  30.5  m  and  inter-element  spacing  of  3.5m.  A  typical  sound  speed  profile  measured  near 
the  trial  site  and  used  in  this  analysis  is  shown  in  Figure  2. 
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Figure  1 .  Arrangement  and  Location  of  Acoustic  Elements  Used  at  the  Experimental  Site. 
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Figure  2.  Sound  Speed  Profile  Measured  Near  the  Trial  Site  and  Used  in  This  Study. 


III.  ANALYSIS  METHODS 


A.         INVERSION  OF  BOTTOM  SOUND  SPEED 

In  determining  the  bottom  sound  speed  a  method  was  performed  to  determine  the 
critical  angle  of  the  seabed  from  the  vertical  directionality  of  the  ambient  noise  in  the  water 
column  (Buckingham  and  Jones,  1987).  The  properties  of  the  ambient  noise  that  we  are 
specifically  concerned  with  are  the  vertical  directionality  and  the  spatial  coherence  of  the 
fluctuations  of  two  sensors  vertically  displaced  in  the  noise  field.  These  two  quantities  are 
related  through  a  Fourier  transform  relationship,  which  is  given  by 

1 
T(Q)  =  -  f  F(6)  exp[iQ  cos(6>)]  sm(0)d0  (1) 

where  G  is  the  polar  angle  measured  down  from  the  zenith  and  F(Q)  is  a  dimensionless 
directional  density  function  which  represents  the  distribution  of  the  noise  power  in  the 
vertical.  F(Q)  can  be  normalized  according  to  the  condition 

i|>(e)sin(e)de  =  i.  (2) 

This  indicates  that  in  the  limit  as  Q  — »  0  the  value  of  the  coherence  function  goes  to  unity. 

Buckingham  and  Jones  provided  a  simple  example  for  isotropic  noise  in  which  F(6) 
is  independent  of  0.  According  to  Eqs.  (1)  and  (2)  the  coherence  function  in  this  case  is 
shown  to  be 


where  Q  is  a  dimensionless  angular  frequency  defined  as 

col 
Q  =  — ,  (4) 

co    is  the  angular  frequency,  cw  is  the  sound  speed  in  the  water  column,  and  /  is  the 
separation  between  the  two  sensor  positions.  This  example  is  shown  in  Figure  3. 


isot 


Figure  3.  Example  of  Ideal  Vertical  Coherence  Function  Versus  Frequency. 

This  function  predicts  that  the  coherence  of  the  noise  at  two  sensors  is  real  due  to 
symmetry  of  the  noise  field  about  the  horizontal.  For  an  asymmetrical  noise  distribution  in 
the  vertical,  the  coherence  function  is  complex  with  a  large  imaginary  component.  A 
symmetrical  noise  distribution  is  also  predicted  by  Buckingham  and  Jones'  theoretical 
model  of  shallow  water  noise,  where  the  coherence  function  for  shallow  water  noise  is  an 
oscillatory  function  that  decays  with  increasing  Q .  At  a  given  frequency,  the  directional 
density  function  of  the  noise  predicted  by  the  model  is  closely  approximated  by 


WH 


a. ,  for 


71 


■0 


<0< 


n 


+  6. 


(5) 


0,  otherwise 
where  a,     is  independent  of  G.     Then,  the  noise  propagates  in  directions  close  to  the 

horizontal,  with  no  acoustic  energy  traveling  with  grazing  angles  greater  than  0C  (the  critical 
grazing  angle).  The  contribution  of  local,  wind-driven  sources  to  the  overall  noise  field  is 
approximately  isotropic,  that  is,  the  directional  density  function  of  the  continuous 
component  of  the  noise  can  be  expressed  as 

F2{Q)  =  a2,  for0<e<7t,  (6) 


where  a2  is  independent  of  0.    By  superposing  Eqs.  (5)  and  (6)  we  have  a  representative 
directional  density  function  for  the  total  noise  field  in  shallow  water: 

/rw(e)=F1(e)+F2(e) 


(a}+a2),  for 


71 


e„  <e< 


71 


+e. 


V- 


(7) 


a2, 


otherwise. 


Figure  4  shows  a  sketch  of  the  vertical  distribution  expected  in  a  shallow  water 
environment.  Through  the  normalization  condition  in  Eq.  (2),  the  critical  grazing  angle  can 
be  expressed  in  terms  of  the  other  two  parameters  as 


sin 


(O-1? 


(8) 


and  the  model  coherence  function  of  the  noise  represented  by  the  distribution  in  Eq.  (7)  is 
found  from  Eq.  (1)  to  be 


Ue)  =  |sin(Q)+^ 


sin[Qsin(0jl. 


(9) 


Using  this  model,  the  spatial  coherence  of  two  vertically  displaced  sensors  of  the 
NW  VLA  was  computed  from  extracted  ambient  noise  data.  By  spatial  coherence  is  meant 


Figure  4.  Model  of  the  Vertical  Distribution  of  Ambient  Noise  in  a  Pekeris  Waveguide. 
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the  time-averaged  crosscorrelation  coefficient  of  the  noise  observed  by  hydrophones 
separated  vertically  in  the  water  column,  i.e., 

PxyM 

meas  ~  run  \  D  i  \ '  (*°) 

where  Pxy(a>)  represents  the  cross  spectral  density  as  a  function  of  frequency  for  two 
separated  phones  and  /^(gd)  and  P^co)  represent  the  spectral  density  for  each  one  of  the 
separated  phones.  To  search  the  two  unknown  parameters  a2  and  0f  in  Eq.  (9),  an 
optimization  procedure  was  implemented  by  finding  a  least  squares  fit  between  the  real  part 
of  the  measured  coherence  function,  Ymeas ,  and  the  theoretical  coherence  function,  Ysw.  By 
defining  the  ambiguity  function 

0  =  ZlrsJQ'fl2?ec)-Real{rmeas(Q)}|2 ,  (11) 

n=l 

values  for  a2  and  0C  are  obtained  when  O  is  a  minimum. 
B.         INVERSION  OF  BOTTOM  DENSITY 

To  determine  the  bottom  density,  the  concept  of  an  effective  depth  for  an  ideal 
waveguide  was  employed  (Chapman  and  Ward,  1989).  By  considering  the  phase  change  of 
a  plane  wave  reflected  at  a  fluid/fluid  boundary,  Weston  introduced  the  concept  of  the 
effective  boundary  depth,  arguing  that  one  could  view  the  reflection  as  taking  place  at  an 
imaginary  pressure  -  release  boundary  located  at  a  specific  depth  below  the  true  boundary. 
He  proposed  that  the  normal  modes  of  the  Pekeris  model  should  be  similar  to  those  of  an 
ideal  pressure  -  release  waveguide  whose  depth  was  equal  to  the  true  water  depth  plus  the 
effective  boundary  depth.  The  normal  mode  wavenumbers  for  the  Pekeris  model  can  then 
be  estimated  easily  without  numerical  iteration. 

The  Pekeris  model  consists  of  a  uniform  water  layer  of  constant  depth  bounded 
above  by  a  pressure  -  release  surface  and  below  by  a  uniform  fluid  half  space  having 
different  density  and  sound  speed.  Consider  a  plane  wave  of  frequency/in  a  fluid  medium 
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of  sound  speed  cw  incident  upon  a  plane  boundary  with  another  acoustic  medium  at  angle  9 
relative  to  grazing,  as  depicted  in  Figure  5. 


Imaginary  Boundary  XLff 

Figure  5.  Geometry  Used  to  Derive  The  Effective  Boundary  Depth. 
The  ray  segments  s  and  s '  are  given  by 


s  = 


and 


A// 
sinG 


2  AH  cos2  9 


(12) 


(13) 


sinG 
and  the  phase  of  the  reflected  plane  wave  at  points  A  and  B  is  given  by 

+A=K/+v  (14) 

and 

$B=2kws+V,  (15) 

CO 

in  which  kw  =  — .  Requiring  §A  =tyB  the  effective  depth  of  the  boundary  is  defined  by 


He(Q)-H 


2A:  sinG 


(16) 
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For  an  isospeed  water  column  with  sound  speed  cw  and  density  p  overlying  a 
homogeneous  solid  medium  with  compressional  sound  speed  c, ,  shear  speed  cs ,  and 
density  p.,  Chapman  and  Ward  (1984),  have  shown  (assuming  cs  <  cwand  bottom 
attenuation  is  negligible)  that  the  effective  depth  is 


H=H  + 


1-2 


'V 


-|2 


\cwJ 


g 


sin  6, 


(17) 


where  H  is  the  water  depth  and  Qc  is  the  grazing  critical  angle.  If  cs  «  cw  then  Eq.  (17)  is 
reduced  to 


H=H  + 


g 


(18) 


kw  sinGf 

The  boundary  conditions  for  this  ideal  waveguide  require  the  pressure  to  vanish  at 
the  upper  and  lower  interfaces.  Correspondingly,  eigenfunctions,  \\i{z) ,  must  satisfy 
(Medwin  and  Clay,  1998) 

v|/(z)|2=o=<|)andv|/(z)L=/fc=0.  (19) 

Solutions  of  the  eigenfunctions  for  these  boundary  conditions  are 


\|/(z)  =  si 


'  wrz  > 


sin 


kHsJ 


(20) 


and  the  vertical  component  of  the  wavenumber,  y,  must  then  satisfy  the  quantization 
condition 


Yn 


h7 


(21) 


where  n  is  the  integer  and  designates  the  mode  number.    The  horizontal  values  of  the 
wavenumber  Kn  are  then  given  by 
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and  the  requirement  that  A^nbe  real  gives  the  mode  cut-off.    By  substituting  y  n  from  Eq. 


(2 1 )  and  kw  =  — ,  Kn  can  be  written  as 
c.. 


K  = 


(     V 
'  CO  ' 


(      \2 
nn 


\c„J 


\HJ 


(23) 


Furthermore,  the  group  velocity  of  the  n*  mode  is  given  by 


dm 


(24) 


therefore, 


2         \(    r.,V  (   „„\ 


u.  = 


CO 


CO 


VCW 


nn 


=  c„,Jl- 


(        ^2 

MIC.. 


(25) 


vco//ey 


The  group  velocity  un  and  the  mode  travel  times  are  related  by 

R 


(26) 


where  R  is  the  horizontal  distance  between  source  and  receiver  and  t„.  is  the  travel  time  of 
the  n"1  mode.  Combining  Eqs.  (25)  in  (26)  the  mode  travel  times  can  be  written 

1 


(27) 


1- 


nc. 


2JHe 


From  the  definition  of  effective  depth,  this  becomes 


tn(fh 


R 


1- 


-i2 


«7isinB, 


^27r/Hsine/ 


+  g 


-iv  J 


-1/2 


(28) 


which  is  the  final  expression  for  the  mode  travel  times  .   Since  we  know  the  real  depth  H, 
and  cw  and  8C  were  determined  from  the  previous  analysis,  then  we  can  invert  for  g  by 
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measuring   the   mode   travel   times   over   a   range   of  frequencies    for   different    SUS 
transmissions.  Assuming  p   «  l,theng«/?,. 

C.         INVERSION  OF  BOTTOM  ATTENUATION 

Mode  attenuation  as  a  function  of  range  can  be  estimated  from  several  computed 
SUS  spectrograms.  The  mode  attenuation  per  km,  an ,  satisfies  the  relation 

ATLn=anAR,  (29) 

where  &TLn  and  AR  are  the  mode  transmission  loss  due  to  attenuation  and  range  between 

SUS  drops,  respectively.  Neglecting  volume  attenuation  in  the  water  column,  the  decrease 

in  mode  amplitude  is  expected  to  be  caused  only  by  bottom  attenuation  and  cylindrical 

spreading  (TLspTeid  =  lOlog/?).  Since  mode  n  can  be  thought  of  as  due  to  rays  propagating  at 

angle  0„ ,  the  mode  attenuation  can  be  considered  the  result  of  multiple  bottom  reflection 

losses. 

The  number  of  bottom  bounces  over  1  km  of  horizontal  propagation  can  be  found 

by  considering  the  angle  of  elevation  or  depression  G  of  the  local  direction  of  propagation 

of  the  waveform.    Referring  to  Figure  6,  the  angle  0  and  the  diagonal  distance  for  one 

bounce  Z,,  are  related  by 

1H 
sin0  = ,  (30) 

A 
where  H  is  the  water  depth.  Defining  the  total  diagonal  distance  L  corresponding  to  R  =  1 
km  of  horizontal  propagation  by 

L  =  ~^z,  (31) 

COS0 

then  L  and  Z,,  can  be  related  in  order  to  define  the  number  of  bottom  bounces  per  km  as 

^=Z=tftan0 
Lx        2H 

The  attenuation  for  mode  n  can  then  be  expected  to  be  due  to  W(0„ )  bottom  bounces  per 
km  producing 

a„  =  M0.K(0.)'  (33) 
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where  0n  is  the  corresponding  propagation  angle  for  mode  n  and  aM(9„)is  the  bottom  loss 
per  bounce  for  a  single  reflection  at  this  angle. 


R=l  km 


Figure  6.  Geometry  Used  to  Find  the  Number  of  Bottom  Bounces. 

The  reflection  coefficient  at  a  fluid-fluid  interface  for  a  plane  wave  traveling  at 
angle  6  is  defined  by 


R  =  pbcbs\nQ  +  ipwb2 
pbcbsw$-ipj>2 


(34) 


where  cw  and  pw  are  the  sound  speed  and  density,  respectively,  in  the  upper  medium,  c,  and 
Pb  are  the  sound  speed  and  density,  respectively,  in  the  lower  medium,  and 


b2  = 


( '     Y 


\cwJ 


-11/2 


cos2e-i 


(35) 
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If  there  is  no  attenuation,  then  |R|  =  1  for  incidence  beyond  the  critical  angle.  However,  if 
the  bottom  is  lossy  then 

|*(9)  =  1-8(0),  (36) 

and  the  bottom  loss  per  bounce  can  be  defined  as 

ajej  =  -201og(j*(ej)=20e(ej.  (37) 

Combining  the  equations  above,  the  reduction  in  the  reflection  coefficient  amplitude  can  be 
determined  by 

v  nJ     2O7V(0„)     20ARN{Qn) 
Bottom  attenuation  can  be  introduced  by  defining  an  imaginary  component  of  the 
bottom  sound  speed, 

cb^cb-ih,  (39) 

such  that  the  wavenumber  in  the  lower  medium  becomes 

kb=  —  -+kb+iae,  (40) 

where 

<*,=—•  (41) 

cb 

Since  the  solution  in  the  lower  medium  goes  like 

eikhR  ^  eikbRe-acR  ^  (42) 

then  ae  is  found  to  produce  the  exponential  decay  in  range.   This  then  causes  transmission 
loss  according  to 

ATLallen  =  8.686cteA/?  =  aAR ,  (43) 

where  a  is  the  bottom  attenuation  in  units  dB/m,  and  can  be  written  as 

a  =  8.686^.  (44) 

cb 

In  order  to  estimate  the  bottom  attenuation  (dB/m),  a  simple  matching  algorithm 
may  be  employed  which  searches  over  values  of  5  (the  complex  bottom  sound  speed)  until 
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the  reduction  in  the  reflection  coefficient  amplitude  agrees  with  the  observed  value  of  c(6n) 
defined  in  Eq.  37.  Once  values  for  8  are  obtained  (for  various  frequencies  and  mode  angles 
0n),  then  Eq.  43  may  be  employed  to  produce  the  final  estimate  for  the  bottom  attenuation 
a. 
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IV.  DATA  ANALYSIS  AND  RESULTS 


A.         INVERSION  OF  BOTTOM  SOUND  SPEED 

As  described  in  Section  III.  A.,  the  inversion  for  bottom  sound  speed  is  based  on  the 
evaluation  of  the  vertical  coherence  of  the  ambient  noise  field.  For  this  analysis,  the 
ambient  noise  was  extracted  from  segments  of  the  SUS  data  sets  between  actual  SUS 
signals  (and  allowing  for  reverberation  to  fade).  The  spectrum  of  the  ambient  noise 
computed  showed  little  variation  from  low  frequency  up  to  600  Hz  indicting  no  strong 
unexpected  noise  sources  were  present.  Also,  the  imaginary  component  of  the  coherence 
between  the  two  hydrophones  was  found  to  be  small  as  required.  This  indicates  that  the 
ambient  noise  field  at  the  trial  site  for  frequencies  below  600  Hz  has  a  symmetrical 
distribution  in  the  vertical  direction  and  the  data  is  suitable  for  inversion  by  the  method 
described  previously. 

It  should  be  noted  that  the  inversion  method  described  in  Section  III.  A.  was 
originally  defined  for  an  isospeed  water  column  overlying  a  homogenous  bottom.  As  is 
obvious  from  the  local  sound  speed  data  (refer  to  Figure  2),  there  is  significant  variation  in 
sound  speed  over  depth.  However,  if  the  sound  speed  change  is  not  too  large  over  the 
vertical  separation  of  the  phones,  the  inversion  method  should  still  be  applicable  by  simply 
adjusting  the  local  angles  to  the  appropriate  values  at  the  bottom  interface  using  Snell's  law. 

Since  the  data  from  the  upper  two  phones  of  the  VLA  was  analyzed,  a  local  sound 
speed  of  cw  «  1485  m/s  was  used  in  the  calculation  of  the  ambiguity  function  (refer  to  Eq. 
11). 

The  fit  between  the  theoretical  coherence  obtained  from  the  model  and  the  real 
component  of  the  observed  coherence  is  shown  in  Figure  7.  By  computing  the  ambiguity 
surface  for  the  optimization  of  two  parameters  (a2  and  0f )  an  optimal  value  was  found  at 

9C * 1 8°  in  the  water  column.    Figure  8  shows  the  ambiguity  surface  for  the  optimization 
with  two  parameters,  a2  and  0C  whose  minimum  corresponds  to  a2 «  0. 1 8  and  mean  critical 
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angle  0f«18°±l°  for  4  different  SUS  explosive  data.  A  grazing  critical  angle  of 
approximately  19°±1°  near  the  bottom  was  determined  from  Snell's  law.  The  sound  speed 
close  to  the  bottom  was  1476  m/s,  and  therefore  the  bottom  sound  speed  c2«1563±10  m/s. 
This  value  was  compared  to  values  obtained  by  a  sophisticated  genetic  algorithm  technique 
employed  at  the  University  of  Rhode  Island  (Miller  and  Potty,  1998).  They  found  values  in 
the  upper  few  meters  of  the  sediment  in  agreement  with  the  result  reported  here. 
B.         INVERSION  OF  BOTTOM  DENSITY 

To  invert  for  bottom  density,  the  method  of  ideal  waveguide  effective  depth  was 
introduced.  However,  this  method  assumed  a  constant  sound  speed  in  the  water  column.  In 
fact,  this  was  fundamental  in  the  definition  of  the  ideal  waveguide  eigenfunctions.  For  that 
reason,  this  method  may  not  be  expected  to  produce  very  good  results.  We  will  attempt  to 
invert  for  density  anyway  by  assuming  the  water  column  is  isospeed  with  nominal  sound 
speed  cw  »  1485  m/s.  Consistent  with  this  value,  we  shall  assume  the  critical  angle  0C  «  18°. 
In  order  to  obtain  separated  mode  travel  times,  spectrograms  were  computed  for  various 
SUS  explosive  data.  The  theoretical  mode  travel  times,  for  the  ideal  waveguide  were  then 
computed  as  a  function  of  range  R  between  source  (SUS)  and  receiver  (VLA),  water  depth 
H,  water  sound  speed  cw  and  critical  angle  Gc.  A  comparison  between  theoretical  and 
measured  travel  times  for  all  1 1  SUS  data  sets  produced  an  average  value  for  the  density 
ratio  of  1 .6.  However,  the  changes  in  the  theoretical  curves  for  various  density  ratios  was 
not  large  whereas  the  mode  arrival  structures  were  fairly  ambiguous.  Thus  the  uncertainty 
is  very  large.  It  is  also  possible  that  the  theoretical  mode  travel  times  have  significant  errors 
due  to  the  mismatch  in  sound  speed  profile.  However,  the  estimated  value  of  pb «  1 .6  g/cm3 
from  the  inversion  method  seems  quite  realistic.  The  results  supplied  by  Miller  and  Potty 
suggest  a  more  reasonable  estimate  for  the  bottom  density  is  about  1.8  g/cm3.  Figures  9a 
and  9b  show  the  mode  travel  times  for  a  phone  from  a  SUS  explosive  data  and  those  for  an 
ideal  waveguide. 
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C.         INVERSION  OF  BOTTOM  ATTENUATION 

A  method  was  developed  in  Section  III.C.  to  invert  for  bottom  attenuation  using 
measures  of  modal  transmission  loss  as  a  function  of  range.  In  order  for  this  method  to 
work  successfully,  the  decrease  in  amplitude  due  to  bottom  attenuation  must  be  observable 
over  the  transmission  loss  due  to  cylindrical  spreading.  For  example,  as  the  range  between 
the  SUS  drop  location  and  the  VLA  increases  from  1 5  km  to  30  km,  the  change  in  TL  due 
to  cylindrical  spreading  is  approximately  3  dB.  This  equates  to  about  0.2  dB/km.  In  order 
for  bottom  attenuation  to  be  distinguishable,  Eq.  (38)  suggests  that  the  quantity  207V  8 
should  be  of  this  same  order  of  magnitude  or  larger. 

A  simple  analysis  assuming  bottom  attenuation  increase  linearly  with  frequency 
shows  that  8  increases  with  increasing  frequency.  Specific  values  can  be  obtained  using 
Eqs.  (34)  -  (44)  with  realistic  values  of  bottom  attenuation  a.  However,  as  frequency 
increases,  the  angle  of  propagation  for  a  specific  mode  decreases  thereby  decreasing  N,  the 
number  of  bottom  interactions  per  km. 

If  one  assumes  a  reasonable  value  of  a  ~  0.05  dB/km/Hz,  a  100  Hz  signal  suffers  a 
0.5  dB/km  loss  due  to  attenuation  for  angles  of  propagation  of  about  12°.  At  this  frequency 
in  this  waveguide,  that  corresponds  roughly  to  Mode  2.  Unfortunately,  this  mode  is  not 
clearly  distinguishable  at  this  frequency  from  the  available  SUS  data.  The  lowest  modes 
can  be  isolated  at  lower  frequencies,  but  then  the  attenuation  is  negligible.  At  higher 
frequencies,  the  propagation  angles  for  these  low  modes  is  near  grazing,  and  again  the 
influence  of  attenuation  is  small. 

The  result  of  this  analysis,  unfortunately,  shows  that  this  data  does  not  avail  itself  to 
inversion  for  bottom  attenuation  using  the  method  developed  previously.  Other  processing 
methods  may  be  able  to  extract  the  necessary  information  in  order  to  utilize  this  method,  but 
will  not  be  explored  further  here. 
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Figure  7.  Coherence  of  the  Measured  Ambient  Noise  between  Two  Hydrophones;  Real 
Component  (Dashed  Curve)  and  Imaginary  Component  (Dotted  Curve).  A  Good  Fit 
Between  the  Theoretical  Coherence  Obtained  from  the  Model  (Solid  Line)  and  the  Real  of 

the  Measured  Coherence  is  Observed. 
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Figure  8.  Ambiguity  Surface  for  the  Optimization  with  Two  Parameters,  a2  and  9C.  The 
Minimum  is  Marked  by  X  which  Corresponds  to  a2  *  0. 1 8  and  0C »  1 8°  for  4  Different  SUS 

Explosive  Data. 
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Figure  9b.  Mode  Travel  Times  for  a  Phone  from  an  Ideal  Waveguide. 
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V.  CONCLUSION 


This  thesis  incorporated  several  methods  which  may  be  applied  directly  to  measured 
data  to  provide  a  relatively  simple  means  of  obtaining  reasonable  estimates  of  ocean  bottom 
parameters  from  minimal  information.  The  analysis  of  the  data  from  the  Mid-Atlantic 
Bight  experiment  recorded  on  the  NW  VLA  produced  the  following  results: 

1 .  Ambient  noise  was  used  to  measure  both  the  vertical  directionality  and  spatial 
coherence  of  the  noise  field.  By  examining  the  noise  data,  we  have  been  able  to 
estimate  Gc  and  hence  the  sound  speed  cb  in  the  sediment.  The  bottom  sound 
speed  deduced  was  approximately  1563110  m/s. 

2.  The  method  of  effective  depth  was  employed  to  estimate  the  bottom  density  pb. 
Although  this  method  was  derived  for  an  isospeed  water  column  environment,  it 
is  expected  to  produce  reasonable  estimates  for  most  shallow  water 
environments.  By  comparing  mode  travel  times  extracted  from  SUS  data 
spectrograms  and  comparing  these  with  values  for  an  ideal  waveguide  with  an 
effective  depth,  the  bottom  density  was  estimated  to  be  pb  »  1.6  g/cm3.  The 
uncertainty  in  this  value  is  very  large,  however,  partly  due  to  the  difficulty  in 
distinguishing  mode  arrivals  clearly  in  the  spectrograms. 

3.  In  order  to  distinguish  bottom  attenuation  from  simple  cylindrical  spreading,  it 
was  determined  that  values  for  mode  transmission  loss  were  required  for  high 
modes  at  high  frequencies.  Unfortunately,  the  SUS  spectrogram  data  could  not 
adequately  provide  this  information,  and  no  estimate  of  bottom  attenuation  was 
produced.  It  may  be  possible  to  use  the  method  outlined  if  a  modal 
decomposition  of  the  data  on  the  VLA  could  be  performed.  However,  due  to  the 
relatively  short  aperture  of  the  array  (compared  to  the  waveguide  depth),  this 
seems  unlikely.  Furthermore,  the  source  level  of  different  SUS  detonations  may 
not  be  repeatable  enough  to  use  this  data  for  accurate  estimations  of  bottom 
attenuation. 
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Comparing  the  inversion  results  for  bottom  sound  speed  and  density  with  values 
obtained  from  genetic  algorithm  techniques  provided  by  the  University  of  Rhode  Island, 
good  agreement  was  found  particularly  for  the  sound  speed.  However,  the  analysis  of  the 
SUS  data  proved  more  difficult  than  originally  anticipated.  Much  of  this  has  to  do  with  the 
uncertainty  in  absolute  drop  location  and  repeatable  source  signals.  Another  part  of  the 
problem  was  the  reliance  on  low  mode  information  which  is  difficult  to  characterize  for 
SUS  data  measured  on  only  eight  hydrophones  with  total  aperture  less  than  half  the  width 
of  the  waveguide.  It  is  likely  that  this  analysis  would  be  more  successful  with  different 
more  well-defined  sources  and  a  full  water  column  array. 
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